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Abstract. We numerically investigate the complex nonlinear dynamics 
for two independently coupled laser systems consisting of (i) mutually 
delay-coupled edge emitting diode lasers and (ii) injection-locked quan- 
tum nanostructures lasers. A comparative study in dependence on the 
dynamical role of a parameter, which determine the phase-amplitude 
coupling of the optical field, in both the cases is probed. The variation 
of a lead to conspicuous changes in the dynamics of both the systems, 
which are characterized and investigated as a function of optical injec- 
tion strength r\ for the fixed coupled-cavity delay time r. Our analysis 
is based on the observation that the cross-correlation and bifurcation 
measures unveil the signature of enhancement of amplitude-death is- 
lands in which the coupled lasers mutually stay in stable phase-locked 
states. In addition, we provide a qualitative understanding of the phys- 
ical mechanisms underlying the observed dynamical behavior and its 
dependence on a. The amplitude death and the existence of multiple 
amplitude death islands could be implemented for applications includ- 
ing diode lasers stabilization. 



1 Introduction 

The semiconductor diode lasers are known to be very sensitive to the external optical 
perturbations such as, optical self-feedback, optoelectronic feedback, optical injection. 
When the diode laser is subjected to optical injection by another diode laser then the 
radiation emitted from mutually delay-coupled diode lasers is well behaving, well 
understandable and well classifiable in terms of complex nonlinear dynamics [2 . On 
the one hand these dynamical instabilities are undesired features and disturb the many 
applications where one needs the constant stable power but one the other hand they 
may allow new methods for the secure communications using chaos synchronization. 
So the systematic study and control of these nonlinear dynamics provide fundamental 
insights into the underlying physics of the system [3]. On the basis of which one can 
redesign the device [TTj or improve the processing, or simply exploit the dynamical 
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performance of a system to ones advantage. Amplitude death [T] is one of the collective 
fascinating phenomena [4], where coupled lasers drive each other to a fixed point 
and stop the oscillatory dynamics. The coupling induced stabilization is known to 
occur either by changing the stability of unstable fixed points which already exists 
in the absence of coupling or the stationary state which can be created newly by the 
coupling. Because of very diverse time scales involved in the dynamics, the routes by 
which the two coupled lasers reach their stable states in the presence of finite-delayed 
interaction remain poorly understood . In this sprit, we have addressed the question 
whether these dynamical behaviors of a coupled diode lasers and particularly quantum 
nanostructures based diode lasers can be controlled via amplitude death and the kind 
of dynamics one encountere during the route to reach the ultimate stable state. 



2 The Model for the delay-coupled diode lasers system 

For numerical simulations, the time evolutions of the complex electric field Ei^ 2 (t) 
and the effective carrier density Ni t z(t) (n t h ~ ^1,2) in the laser medium for each 
delay-coupled diode lasers are modeled. 

^ = (l+ia)N 1 (t)E 1 (t) + r,e- iU2T E 2 {t-T), (1) 

T^ = Ji-iVr-^ + lH^WI 2 , (2) 

= (1 + ia)N 2 {t)E 2 {t) + V e-^E^t - r), (3) 

T^ = J 2 -N 2 -[2N 2 + l]\E 2 (t)\\ (4) 

In order to analyze the complex dynamical behaviour of two mutually delay-coupled 
diode lasers system (as shown in figure 1), the LK equations can be written in a 
standard normalized form [2] : where rj is the effective coupling strength, i.e., the 
fraction of light of one laser injected into the other laser and vice versa, Ji^'s are the 
effective current densities (with the threshold value subtracted out), T is the ratio of 
the carrier lifetime to the photon lifetime, a [7] is the linewidth enhancement factor 
or phase-amplitude coupling factor, and r is the coupled cavity time taken by the 
light to cover the distance between the two lasers system. u>i t2 are the optical angular 
frequencies of the free running diode lasers 1 and 2. Ei, 2 (t — r) are the fields delayed by 
one coupling time r = L/c and uj\^ 2 t are the coupling phase detuning. For simplicity, 
we have taken two similar diode lasers for which the frequency detuning between them 
is assumed to be very small or nearly zero. In order to mimick the real experimental 
situation, we have introduced some small noise in our model. We have done numerical 
integration of the above equations by using Runge-Kutta fourth-order method with a 
step size = r/n, where n — 1000 is chosen based on the response time of photodiode 
in experimental setup[2 ]. The dimensionless parameters are taken as J\ :2 = 0.165, 175 
and T = 1000 [2j. In order to scan and understand the different dynamical regimes, 
we use equation 5 for correlation measurement pQ .We observe the route to amplitude 
death of low frequency complex dynamics in the output power of a diode laser (slave) 
when subjected to optical injection from another diode laser (master). In the context 
of a system of two delay-coupled diode lasers and Quantum dot nanostructures lasers, 
we emphasize the effect of a on the control of complex dynamical instabilities near 
the phase flip transition regimes as a function of coupled-cavity time delay r and the 
optical injection strength r\. Shifting of different phase-correlated dynamics, such as 
phase- flip bifurcation [2] and strange bifurcation [1] was observed in our previous work 
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Fig. 1. (Color online) Schematic diagram of a Delay-coupled diode laser system. 



when T) is varied for a particular r. In this present work, we show that these phase flip 
transitions does not occur abruptly at a particular value of rj and a. Instead we find 
a coupling strength region around the phase flip transition where the co-existence of 
multi-attractors occurs as shown in figure 2 and figure 3. We show that the phase 
flip bifurcation occurs from in-phase amplitude death to anti-phase amplitude death 
and the strange bifurcation occur from anti-phase to in-phase transition regimes |10j . 
The existence of multiple attractors near the regime of strange bifurcation (as shown 
in fig. 2. and fig. 4.) has raised the issue of whether the a plays a crucial role or not. 
We also study the effect of r\ along with the phase-amplitude coupling factor on the 
dynamics in the amplitude death regime and extends the study for nanostructures 
quantum dot lasers. One of the Key features of semi-conducteur lasers is their 3 — dB 
modulation bandwidth, which is limited due to the presence of strong relaxation 
oscillation damping rates. This work also theoretically focuses on the impact of strong 
injection in a quantum nanostructure semiconductor laser with large variation of the 
phase-amplitude coupling factor through the non-linear dynamics and the modulation 
response at zero-detuning. The combination of the strong injection, optimized phase- 
amplitude coupling factor and the zero-detuning case shows the possibility to reach 
the stable state, which is fully suitable for laser control applications. 



3 Influence of a in delay-coupled edge-emitting diode lasers and 
injection-locked quantum nanostructures lasers 

The linewidth enhancement factor a describes the coupling between the real and 
imaginary parts of the material refrective index of the gain medium, and is given by 
the ratio of their derivatives with respect to the carrier density. Its value is known 
to depend on the carrier concentration, photon wavelength, inhection current and 
temperature. For quantum nanostructures based diode lasers (quantum dot lasers ), 
the value of a can in fact be controlled by the desirable design of the device structure 
or by the filtered optical feedback or injection techniques. Unlike the quantum dot 
lasers, the a in bulk semiconductor lasers can be modified just by varying the specific 
combination of control parameters in a particular range. In order to analyze the role 
of a in the various dynamical regimes of mutually delay-coupled diode lasers, we use 
a normalized cross-correlation function defined as 



Note that the measured output powers of the lasers, Pi ,2 = |^i,2| 2 do not explicitly 
depend on the optical phases of the electric fields. We wish to probe the role 
of phases of optical fields through a in the intensity cross-correlation near the phase 



C = 



((P 1 (t)-(P 1 (t)))(P 2 (t)-(P 2 (t))) 



(5) 



V((Pl(t)-(P 1 (tW)((P2(t)-(P2(tW)' 
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Fig. 2. Plot of the cross-correlation C(t) versus a for a fixed time delay r = 14 (in units of 
cavity photon lifetime). 




bifurcation or transitions of two mutually delay-coupled diode lasers. We start by 
analyzing the modulation properties of the coupled laser system from the rate equa- 
tions (1 — 4). In Figure 4, the symbol AD represents the signature of ultimate death 
state while BS and MS represent the Bistability and multistability between in-phase 
amplitude-death island and anti-phase amplitude-death island respectively. Recall 
that the phase-space plot of this ultimate death state is represented as AD in hgure 
5, and the MD dynamics is shown in figure 6. From the cross-correlation measures, we 
have observed a unusual phase transitions from anti-phase to inphase state (as shown 
in fig. 3) via strange bifurcation as the control parameter, a is varied for a particular 
value of r\. So the phase-amplitude coupling factor induced multistability has been 
observed in the system as shown in fig. 4. The amplitude death attractors is shown in 
fig. 5 within the AD regimes of fig. 4. The three different attractors (MS regimes in 
fig. 4.), depending on the initial conditions, are shown in fig. 6. as the rj scaned back 
and forth. Injection-locking of semiconductor lasers is one of the most attractive 
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Fig. 4. Bifurcation Plots of laser output powers Pi (open circles) and Pi (filled circles) 
versus coupling strength r\ for a fixed time delay r = 14 (in units of cavity photon lifetime). 




Fig. 5. Plot of Amplitude death attractor of the AD state shown in figure 4 in the parameter 
space of P(t), N(t), 4>(t) for a fixed time delay r = 14. 



research topics since this method induces superior improvement in the high-speed 
characteristics of directly modulated lasers such as increasing the modulation band- 
width, suppressing nonlinear distortion, relative intensity noise [5J f mode hopping 
and reducing chirp [5] . Previous work has focused on realizing high modulation band- 
widths and associated design strategies, analyzed the modulation properties of the 
coupled system in the spectral domain, and numerically investigated the modulation 
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Fig. 6. (Color online) phase space plots of the multi-attractors dynamics (MS state in figure 
4.) between two stable locking states as the initial condition are varied for fixed -q, r and 
initial transients of 10 5 data points are discarded. 



response of the injection- locked system [3]. In order to understand the limiting factors 
in an injection-locked system it is important to investigate the governing theory that 
can be obtained by properly modeling the impact of characteristic parameters such as 
the so-called linewidth enhancement factor ( ev-parameter) . A dimensionless volume 
averaged normalized approach to theoretically evaluate the nonlinear dynamics as a 
function of the injected field ratio and/or the detuning frequency for varied slave laser 
bias cases can be described as follows [8]: 

dY 

— = ZY + eY(Y 2 -P) + r; cos 6, (6) 
dr 

!=L=aZ- ae(Y 2 - P) - ± sin 9 - AM, (7) 
dZ 

T — = P - Z -Y 2 (1 + 2Z -2eY 2 + 2eP) + 77 cos 6>, (8) 
dr 

with Y, describing the normalized field magnitude and Z the normalized carrier den- 
sity. The T-parameter is the ratio of the cavity decay rate to the spontaneous carrier 
relaxation rate. Parameter P is proportional to the pumping current above threshold 
while coefficient e accounts for the nonlinear carrier contribution to relaxation rate. 
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Fig. 7. Calculated bifurcations diagrams of a quantum dot laser with a variable linewidth 
enhancement factor. 



The detuning and phase offset between the master and the slave are denoted by AQ 
and 9, respectively. The normalized injection strength is rj — r]o(P/'y c ) 1 ^ 2 'with 7 C the 
cavity decay rate and 770 the maximum injection strength. In solving the coupled 
normalized differential equations, the normalized field magnitude Y is not at steady 
state, and is thus represented as a dependent term in the normalized field magnitude 
and phase rate equations. In what follows, this model is used for the stability analysis 
of a quantum dot laser. The quantum dot laser under study is a ridge waveguide 
with 500-m cleaved cavity length. The ground state (GS) emission wavelength is at 
1560nm 6 .Because of the carrier filling in the lasing and non-lasing higher energy 
levels, gain compression effects are much larger as compared to bulk or quantum well 
materials [5]. Consequently, it has been shown that the a- factor in quantum dot lasers 
cannot be considered as a constant parameter since it strongly varies from one laser to 
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another and also with optical output power. Figure 7 shows the bifurcation diagrams 
calculated at a constant pump current but for different values of the a-parameter 
ranging from 1 to 15(assumingZ\J7 = 0). The objective of the calculations is to show 
the effects of a large a-factor on the lasers stability. Numerical results point out that 
taking into account such variations reveals strong modifications in the bifurcation 
diagram. On one hand, at low a (case (a)), the laser is always stable while for a = 3 
(case (b)), period one oscillation starts occurring. On the other hand, for a > 3 (cases 
(c), (d), (e) and (f)),the bifurcation diagram exhibits chaos (at low injection ratio) 
followed by a cascade of periodic regimes converging to a stability area (at very high 
injection ratio). From a general point of view, simulations reveal that the smallest the 
a-parameter, the better. As far as the engineering point of view is concerned, these 
results are very important for the optimization of the microwave properties. Although 
the optical injection is used to purify the relaxation frequency as well as the modula- 
tion bandwidth, these numerical results point out that a large a-factor is detrimental 
for the lasers stability and can induce severe degradations in the microwave proper- 
ties. As a conclusion, in order to maintain a wide stability area with optical injection 
associated to good microwave properties, a low a-factor is mandatory in quantum 
dot laser. 



4 Conclusion 

This new concept stabilizes the laser emission by directly controlling the nonlinear- 
ity and, consequently, enhances the stability properties of the system. In general, 
semiconductor lasers with a sufficiently low a would be most interesting for prac- 
tical applications due to the possibility of chirpless operation, and the insensitivity 
to delayed optical feedback or injection. The different dynamics and the strange bi- 
furcation among them is investigated as a function of coupled-cavity time delay r 
and the optical injection strength rj. The correlation measure gives the signature of 
variation in amplitude death islands of complex dynamics of the delay-coupled lasers. 
The shrinkage of in-phase death state and enlargement of out-of-phase death state 
are observed (as shown in fig. 2 and fig. 3) and analyzed when a and r] are varied. We 
provide detailed information about the effect of the variation of a on the dynamics 
of the system over wide ranges of experiment relevant parameters. In particular, we 
give numerical evidence that the stability of the system increases with decreasing a. 
This last point is particularly predominant in a quantum dot laser for which a larger 
a-factor is usually observed. 
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